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Abstract 

We present in this work a generalization of the solution of Gorenstein and Yang for a 
consistent thermodynamics for systems with a temperature dependent Hamiltonian. We 
show that there is a large class of solutions, work out three particular ones, and discuss 
their physical relevance. We apply the particular solutions for an ideal gas of quasi-gluons, 
and compare the calculation to lattice and perturbative QCD results. 
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1 Introduction 



Lattice QCD suggests PQ El E] that at sufficiently high temperature T and/or quark chemical 
potential fi, the strongly interacting matter exhibits a transition from a hadronic phase to a new 
state, where matter is described in terms of fundamental gluon and quark degrees of freedom 
[1], the Quark-Gluon Plasma (QGP). These lattice simulations suggest that the transition 
occurs for temperatures around ~ 190 MeV [3j. Heavy-ion collisions at RHIC, and in the 
future at the LHC, provide us essential experimental data to the quest of the QGP at high 
temperatures and small chemical potential. This experimental search of the QGP needs reliable 
theoretical estimates of various quantities, such as pressure, entropy, deconfinement temperature 
and equations of state (EoS). 

Perturbative QCD at finite T and fi is one of the theoretical tools to compute the various 
QGP quantities. However, strict perturbation theory, which have been pushed up to gg\n{l/gs) 
[5], is reasonable only for extremely high temperatures: at temperatures near Tc it is, in principle, 
not applicable, and further treatments appear to be necessary. Moreover, the perturbative series 
seems to be weakly convergent [5l[6]. Specifically, it is expected that when T ^ Tc, the plasma 
behaves like an ideal gas of quarks and gluons, but the perturbative series goes very slowly 
towards this expectation. A way to attack this problem is through the reorganization of the 
perturbative series. For instance, there are attempts using resummation based on the Hard- 
Thermal-Loop (HTL) effective action [7], in the form of so-called HTL perturbation theory 
[8], [9], and also attempts based on the 2-loop ^-derivable approximation [TOl [11]. The latter 
approach, which assumes a massive quasi-particle formalism [121 [13], leads to results which 
agree remarkably well with lattice data. 

The quasi-particle model of the Quark-Gluon plasma (qQGP) is a phenomenological model 
which assumes non-interacting massive quasi-particles which, with the aid of few parameters, 
like thermal masses, is able to fit lattice QCD data over a wide range of temperatures [I ^[T5l[T6] : 
not only at extremely high temperatures as in strict perturbation theory, or when T > 20Tc 
as in the HTL perturbation theory, but also near Tc. Quasi-particle models intend to describe 
deconfined matter from Tc up to T ^ oo. In the usual quasi-particle model for the quark-gluon 
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plasma fTH [Ti] , statistical mechanics for an ideal massive gas is used, where the mass of each 
quasi-particle is dependent on the temperature (and on the chemical potential) of the gas. 
However, when this idea is naively applied there appear inconsistencies in the thermodynamics 
of the system [T7j. The in depth study of a general solution for the thermodynamics of a system 
composed of particles with thermal masses is the main objective of the present work. 

In Sec. 2 we will formulate the requirements for thermodynamics self-consistency through the 
use of statistical mechanics in the canonical ensemble, and explicitly compute three solutions. 
In Sec. 3 we will use one of these solutions to solve the problem for an ideal gas of quasi-gluons, 
and obtain the generalized thermodynamics relations as a function of the temperature. In Sec. 
4 we discuss which solutions have physical meaning, and in Sec. 5 we analyze and summarize 
our results. 

2 Statistical Mechanics of Canonical Ensemble with a 
Temperature Dependent Hamiltonian 

The partition function of a classical system in a canonical ensemble whose Hamiltonian depends 
on an extra parameter T, which will be identified with the temperature of the system, has 
the same formal structure of the partition function of a system that has a T independent 
Hamiltonian. Specifically: 



where (3 = T~^, is the number of particles, H is the system Hamiltonian which depends on 
the momentum p, the coordinates q of each particle, and the temperature of the system. The 
integral is computed over all momenta and coordinates of the particle (we use the notation 



where A is a thermodynamic function to be determined later. For instance, in the Standard Sta- 





Jdpdp = Jd^^pd^^q). 

The thermodynamics of the system is obtained from the partition function 



A{V,T) = -T\nQr,, 



(2) 
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tistical Mechanics (SM), one works with a temperature independent Hamiltonian and A{V,T) 
is the free energy. 

In the SM, all the other thermodynamics functions can be found from A[V,T) using the 
thermodynamics relations: 



• ■ - (S), 

U = A + TS, (3) 



T 



V 



where P is the pressure, S is the entropy, and U is the internal energy. 

Returning to the case of a temperature dependent Hamiltonian, it can be deduced from 
Eqs. ([1]) and ([2]) the following identity 



7V!(27r)= 

Differentiating it with respect to (3 on both sides, and averaging the resulting expression, one 
obtains 

AiV,T^-T?^-iHi,,,,T)HT{?^^)^0. (4) 

The symbol {g) means the average of the function g in the SM, and it is defined as 



(^7) = (27r)3^iV! J dqdpg{q,p)p, 
where p is the ensemble distribution function, given by 



P 



Qn 

In the SM the last term of the Eq. (jl]) vanishes, and the internal energy is assumed to 
he U = (H). Then, comparing equations ([3]) and (jlj), one concludes that A has to be the 
free energy, as stated before. As the last term of Eq. (jl]) does not vanish in general, the usual 



definition can not be used here. Thus, some questions appear for a temperature dependent 
Hamiltonian: which thermodynamic function is A7 Has U the same definition as in the SM7 

We will show that there are innumerable possible answers to these questions, with the 
constraint that one always has to recover the standard Statistical Mechanics in the limit of 
a T independent Hamiltonian. For any of the possible solutions, it is necessary to redefine 
the connection between statistical mechanics and thermodynamics, in such a way that the 
thermodynamics of the system is also consistently built. 

As Eq. dl]) is the one where the problem is explicit, one should start from it. First, one can 
redefine the free energy as a function of the SM free energy A: 

A'{V, T, /(T)) = A{V, T, /(T)) + aB{V, T), (5) 

where B = B(y, T) is an extra term which will be chosen to let the thermodynamics formulation 
consistent, a is an arbitrary real constant, and /(T) is the additional temperature contribution 
to the thermodynamics functions, which appears because of the T dependent Hamiltonian. In 
the limit of a T independent Hamiltonian, B{T) has to be zero, so one recovers all expressions of 
the SM. The redefined free energy A' has to satisfy all the thermodynamics relations, meaning 
that when replacing Eq. ([5]) into Eq ([3]), one finds: 

^ ■ -(¥),.■ 

and 

A' = U -TS + aB. (8) 

The quantities without prime are those calculated in the SM, i.e. they are calculated regarding 
/(T) as a constant. In order to let Eq. ([8]) to have the same form of the third thermodynamic 
relation given by Eq. ([3]), it is necessary to redefine U and S: 
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(6) 
(7) 



S'{V, T, /(T)) = S{V, T, /(T)) - 7^^^, 

[/'(V, T, /(T)) = T, /(T)) + r?i?(\/, T), (9) 

where 7 and are arbitrary constants, but constrained by a = 7 + 77. Using the second relation 
of Eq. ([3]) and Eqs. ([5]) and ([9]), one obtains 

B _ dB df fdA\ 
and manipulating the partial differential equation, one finds two possible solutions. 



T df (dA\ , , 

B = -i^hrr " = 0. (10) 



idT \df 



V,T 



Using the definitions of U\ S' and A' in Eq. (jl]), one obtains 



dB ^IdH' 

or, 



B = -I — ) a = 0. (Ill 
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\9T 



As the entropy in the SM can be written as 

S = (3{H) +\nQi, = -(Inp), 

so, from Eq. ([9]), one is able to find the connection between the ensemble distribution function 
and the redefined entropy: 
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S'iV,TJiT)) = -(Inp) -7^^. (12) 

With the redefined thermodynamics functions, one then constructs a general model to treat 
the case of a T dependent Hamiltonian, where all the thermodynamics relations are satisfied. 
Hamiltonian functions which depends on the temperature of the system appears in problems 
with mean field approximation, as in the theory of nuclear matter [18j, or in phenomenological 
models of Quark Gluon Plasma |T2l[Ill[IZ|. One is able to recover the thermodynamics consis- 
tency of systems with a T-dependent Hamiltonian if proper care of the extra term B is taken. 
The natural question is: what is the meaning of B7 As the A^! factor in Eq. can not be 
justified classically, the meaning of B is to be interpreted with the aid of quantum mechanics. 

For a quantum system, the canonical ensemble with a Hamiltonian operator H has the 
following definition: 



{g) = Trigp), (13) 

where g is a. given hermitian operator and p is the ensemble density operator. Based in the 
redefined classical internal energy Eq. (Q, one is able to write the redefined quantum internal 
energy, with the help of Eq. (JT3D , as 

U' = {H) +r]B = Tr[p{H + r]B)], 

with B = Bl, where 1 is the unitary matrix. Comparing this equation for U' with Eq. f[T3|) . 
one is lead to assume the term inside the brackets as the total Hamiltonian 



Ht = H + Eo, (14) 

where Eq = 7]B. As B does not depend on the momenta and coordinates, the density operator 
for a general Hamiltonian Ht has exactly the same form as that given by Eq. 0131) . Also, as 



p does not change for the redefined Hamiltonian, the entropy does not change as well. The 
quantum statistical mechanics relations are then written as 



S' 



-(Inp) -7— 
-T\nQl + -iB 



A' 



U' 



(Ht). 



(15) 



Here, Qjj is the partition function written in terms of Ht- Note that for 7 = 0, the ther- 
modynamics relations have the same form as those of the SM. The interpretation for B was 
first given by Gorenstein and Yang in Ref. |l7j, where they make the observation that in the 
standard case of a T independent Hamiltonian, the zero point energy is a constant and it is usu- 
ally subtracted out because experiments measure only energy differences. In the quasi-particle 
model, the dispersion relation is T dependent, and so is the zero point energy of the system. 
Thus, it can not be discarded from the energy spectrum. In this sense, rjB is the system energy 
in the absence of quasi-particle excitations, i.e. the system lowest state energy. Note from Eqs. 
([TU]) . ffTTl) and ffT51) . that when one of the three constants, a, •y or rj are zero, the corresponding 
thermodynamics relations are independent of the two remaining constants, i.e., they cancel each 
other. 

In this section, it was determined all possible mathematical solutions for the formulation 
of a consistent thermodynamics for systems with a T dependent Hamiltonian. The constants 77 
and 7 can assume any value, consequently a as well, but some of the values of these constants 
are of practical use in physics. In the next subsections, three of these particular situations will 
be developed. 



The first solution to be dealt with is the one with 7 = and a = rj, which implies that entropy 
is not changed, while the free energy and the internal energy are: 



2.1 Solution 1 
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U' = -^Tr (Htc-^^'A = -^Tr (i^e"^^) + B. (16) 
Qn Qn ^ ' 

From Eqs. ffTOj) and f|TT]) . one determines the B term: 

where i?o is an integration constant. Using this equation for 5, the redefined thermodynamics 
functions are 




This solution was first used by Gorenstein and Yang [17J to study a gluon plasma with 
a T dependent Hamiltonian, where the gluon dispersion relation was assumed to have a T 
dependent gluon mass. In that case, the /(T) function is given by the square of the quasi-gluon 
mass. This solution was used in several works, where the quark-gluon plasma is treated by a 
quasi-particle model [121 [El [13 [S] , and it can be physically supported: if one views a hadron 
through a bag model, then one deals with a "bag pressure" and a "bag energy" . When heat is 
added to the system, there will be some reminiscence from the bag on the energy and on the 
pressure: B plays this role - generally B is referred to as bag energy or bag pressure. 

2.2 Solution 2 

One other possible solution is = and = 7. Then, according to Eq. (ITSll . the internal energy 
is unchanged, and the entropy and the free energy are given by: 



A' 
5" 



-T In Tr (e"^^^) + = -TlnTr (e"^^) + B 

-^Tr (e-f^^^ In pt) - (3B = -^Tr (e-^^lnp) - pB. 
Qn Qn ^ ' 



(19) 



Eqs. f lip I) and (fTTj) can be rewritten for rj = as: 



dTKTj TdT\dfjy^ T\dT /' ^ ' 

Computing B from Eq. fl20l) . the redefined thermodynamics functions are 

V = {H). (21) 

This second solution was used in Ref. [IS]. It describes the quark-gluon plasma by a quasi- 
particle model satisfying all thermodynamics relations as well. However, here the physics moti- 
vation is that the whole interaction energy goes to the quasi-particle mass, i.e. there is no extra 
term in the usual internal energy expression. 

2.3 Solution 3 

A third solution consists in a = and t] = —7. This is simpler than the others, since with 
a = 0, Eqs. flTUl) and flTTl) for B are not partial differential equations. In this solution, the free 
energy is unchanged, while the entropy and the internal energy are: 

S' = -^Trle-f^^^ In pT)+/3B=^Tr(e-^^\np)+pB, 
Qn Qn ^ ' 

U' = -^Tr (Hre-f'^A = -^Tr (i^e"^^) + B. (22) 
Qn Qn ^ ' 

The redefined thermodynamics functions are then 



s' - -an.) -(f). 

A' = -TlnQjv, 

- - m- 
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Note that the internal energy in Eq. fl23l) is the only one of the three solutions in which the 
connection between the partition function and U' has the same form as in the SM case j20j : 



In the quasi-particle model context of the quark-gluon plasma, the third solution can be 
interpreted as the following: the quasi-particle changes the zero point energy, such that it 
becomes a function of T (just as the quasi-particle mass). This effect occurs because the whole 
interaction energy can not be accommodated in the quasi-particle mass, and part of it is used 
by the vacuum to modify its zero point energy, represented here by B(T). In this solution, 
the entropy definition is also modified, which may be seen as a problem. However, as it is well 
known, the entropy of a system is determined up to a constant that is usually subtracted from 
the system entropy. This procedure can not be done here, because for the case of a T dependent 
dispersion relation, the additive term to the entropy will not be a constant, but also a function 
of the temperature. 

3 Ideal Quasi- Particle Gas Model 

In the previous section, we have constructed a self-consistent thermodynamics for a system 
with a T dependent Hamiltonian. In this section we will apply this theory for one specific case: 
a plasma composed only by gluons, with a vanishing chemical potential. The quasi-particle 
model implies that we will treat the plasma as a non-interacting gas of massive, temperature 
dependent, gluons. This problem has already been solved for two particulars solutions, solution 
1 [T7j and 2 [19] of the previous section, of the general solution given there. To exemplify our 
method we choose the easiest solution, solution 3, which gives the simplest extra term B. We 
will see that we can obtain an algebraic solution for the thermodynamics function, differently 
from solutions 1 and 2. 
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3.1 The a = Solution 

It is interesting to work in a grand canonical ensemble because the thermodynamics functions 
are easier to compute there, but we have developed the theory in a canonical ensemble. Never- 
theless, for the II — case it is not necessary to redefine the whole theory in the grand canonical 
ensemble, because all that is required for this change of ensemble is to use the grand canonical 
partition function, Z, instead of Q^. This connection is justified because both ensembles are 
equivalent in the thcrmodynamical limit, and for /i = the additional thermodynamics relation 
for N, N — —dA/d/i, disappears. The grand partition function is given by: 



Z{z,V,T)= J2z''QNiV,T), (24) 

N=Q 



where we introduced the fugacity z, 



z = e^". (25) 



The distribution function p is 



Z ' 



and the mean value is defined as 



(A) 



E^=o^^^Qiv(V^,T) 



j:N=o^''QNiv,T) ■ 

Using the results given in the Sec. 2.3, one has the redefined thermodynamics functions: 



P'V = T\nZ. 

where the thermodynamic relation A — —PV was used. Prom quantum statistical mechanics, 
the partition function for a free gas is 
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k 



or 



v_ 

where the dispersion relation is uj"^ = k"^ + m?{T). Therefore P' and U' are given by Eqs. fl26|) . 
and take the form 



where u is the gluon degeneracy factor, e'(T) = ^ is the energy density and f{k) = {e^'^ — 1)^^ 
is the Bose-Einstein distribution function. Note the second term in the expression for the energy. 
It exists only because m = m(T), and it will be denoted by b(T) = B(T, V)/V. 
Using the relation S'T = U' + P'V, one obtains: 



dkf{k) 
dkfik) 



_k^ 

F + m2(T) 



k'^ + m?{T) 



(27) 
(28) 



s'iT) = —— / dkf(k)k^^ ^ ' -—— / dkf(k)^= 29 



If one uses the thermodynamic relation S" = —[d{P'V)/dT]v to calculate the entropy 
density, with the help of P' given by Eq. (127|) . one finds exactly the same expression Eq. (12^ 
for s', confirming that the calculation is self-consistent; the additional term b{T) in e'(T) and 
s'(T) guarantees the thermodynamical consistency of the qQGP expressions. If one compares 
this extra term with the usual extra term of qQGP, the Gorenstein-Yang solution (Eq.(5) of 
Ref. [H]), one finds the relation b{T) = dB{T)/dT. This implies that it is necessary to compute 
an integral in T in order to get the bag constant B{T). In the present approach, there is not an 
integral in T, what facilitates analytical calculations. Also, as m(T) depends on the coupling 
constant, in the form of a logarithm in T, B{T) is harder to compute. 
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3.1.1 Pressure 

The pressure P'{T) is given by Eq. fl27|) . and it can be rewritten as 
where r = ^ and Is is Eq. (1551) of Appendix A. Eq. (130|1 becomes 



P'(T) 



TT 



4 2 4 

T^m^ + —Tm^ + — 

90 24 12 25 



^°^4;^ + ^^-4 



1 - m2("+2) 

^ n=l 



with a„ 

gas pressure, Pq 



(-l)"(2ra-l)!!C(2n+l) 
(n+2)!23"+l7r2" 



90 



. Rewriting the above integral as a function of the ideal massless 
one obtains the pressure for the massive gas of quasi-gluons: 



P'(T) = Po 



15 /m\ ^ 15 / ^ 
1- — (- + ' 



47r2 \Tj 27r3 VT 



45 f {, ^ 

le^lrj r^4^ + ^^-4' + 



45 ^ /m\2(«+2)' 



47r4 



n=l 



T 



(31) 



Note that the functional form of m(T) is unknown, i.e. statistical mechanics does not provide 
the function m(T). This problem will be treated on the next section. 

3.1.2 Energy 

The energy density e'(T) is given by Eq. fl28l) . and using the same technique as used in the 
calculation of the pressure, one is able to get the expression for the energy density dependent 
on the temperature: 



e (^) = 7TZ2 / dx^===—j^=^— + ^\m - -—- I dx- 



27r2 7o + r2 e^^^^ - 1 27r2 \^ 2 dT J Jo y/x^ + e^^^ - 1 

where /s and Is are Eqs. and (155]) of Appendix A. Solving these integrals one obtains: 
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f 5 /m\2 15 /m\^ / m 1\ 15 ^ , /m\2("+2) 



5 



3 m 



3 



m 1 



3 ~ ^ ^ /m\2(«+i) 
2^ a„(n + 2) 



n=l 



.(33) 



where Cq = 



30 



is the energy density of the ideal gas. 



3.1.3 Entropy 

As the expression for the pressure for the massive gas of quasi-gluons has been obtained, it is 
easier to compute the entropy density through the use of a Maxwell relation, and one gets: 



s'{T) = 



dP' UTT^ d 



dT 90 or 



15 /m\^ 15 / m\ 



— 



47r2 \TJ \T J 



Using 



45 fmV 



IGtt^ 



45 



m 3 



2(n+2)^ 



47r4 



n=l 



T 



d / m\ dm 1 m 



1 dm^ 



m 



dT\T) dTT r2 2mT dT T^' 



s'{T) is written as: 



15 / 
8^ [t. 



45 



647r4 \T 



m 



45 



n=l 



2(n+2) 



15 



1 - 



3 m 3 



TT 



T 27r2 r2 



\ m 1\ 3 ^ , /m\ 



T 



(34) 



where sq = ^^^i^ ^l^® ideal gas' entropy density. 



3.1.4 Results 

The temperature dependence of the thermodynamics functions of a quasi-particle gas can not be 

completely described using statistical mechanics only. In particular, to determine the expression 
for m(T) it is necessary some further information. First, we take HTL perturbation theory 



15 



[TUt [TB] as the additional information. HTL provides a relation between m{T) and T that, 
asymptotically, is 



m2(T)=iV,^, (35) 

where gs is the strong coupling constant. With this input, one can compute the explicit temper- 
ature dependence of the thermodynamics functions. As we are interested not just in the leading 
order behavior but also in the higher order corrections, we will implement these higher order 
corrections to the thermal mass of the HTL approach [TU], 

Sm'' = -N.^Tmo, (36) 
where uid = J^^dsT is the Debye mass. For the values of the strong coupling constant 



at high temperatures, gs ^ I, expression fl36|) should describe correctly the next-to- leading 
order temperature dependence of the mass. However this equation runs into problems, giving 
a tachyon mass for the gluons, depending on the value of gs- As the aim is to construct a self- 
consistent model to describe the deconfined QCD phase using quasi-particles, we can not take 
this next-to-leading mass as the complete quasi-gluon asymptotic mass but as an approximation 
to the complete mass. Analyzing this expression inside the HTL picture [10], it is seen that 
Eq. (l36ll comes exclusively from hard momenta corrections, i.e. soft momenta corrections do 
not enter. The complete description of the asymptotic mass term requires both soft and hard 
corrections, and these will be simulated through a quadratic gap equation given by, 

m\T) = N^^-Y - ^glTm{T). (37) 

Using Eq. (!37|) in Eq. (13T1) . one is able to compare the present results with lattice QCD data 
[2]. For the calculation of the temperature dependence of the coupling constant, we used [22] : 

2 487r^ 

^ llAr,log[Ar/T, + r,/T,]2' 

with A and Tg phenomenological parameters, and Ai'c = 3 is the number of colors that will be 

used hereafter. The optimal fit is achieved using A = 5.182, T^/Tc = —0.197, and the degeneracy 

16 



factor V = 2{N^ — 1) = 16. In Fig. [T]is plotted the thermodynamics functions Eqs. (PTi) . ([2SD 
and flMl) . using the mass relation Eq. fl37|) . The calculated curves describe quite well the lattice 
data [2]. For comparison, in Ref. [14J, where the qQGP was based on the Gorenstein and Yang 
solution, it was necessary 4 parameters to get an optimal fit of the lattice data: A, Tg/Tc, the 
integration constant Bq and the degeneracy u. 
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Figure 1: Plots of the pressure, energy density and entropy density in our model, and the 
extrapolated lattice data [2j for the pressure (star), energy density (triangle) and entropy density 
(diamond), as a function of T/T^. 



3.2 Asymptotic Behavior of the a Solutions 

In the last section, the a = solution was studied in some detail using explicit expressions. 
For the a ^ cases, however, the calculation of the B term involves non trivial integrals. 
Nevertheless, if one is interested only in the asymptotic behavior of the gluon plasma, an 
explicit calculation can also be done. In this sub-section, we will focus on this calculation. 

The general expression for the a 7^ solution for B is given by Eq. (fTTI) . Using the quasi- 
gluon dispersion relation, one can rewrite B as 



B = BqT^ — - / dTT-^ — / dk , ' =. 39 

aAn^ Jtc dr Jo Jk'^ + m?{T) 
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To analyze the high T behavior of B{T), we manipulate Eq. fl39l) in order to use the form found 
in Appendix A: 



B 



BoT-- 
B,Ti- 





f 


drr 


xdim?{T) 


r(3)r^/, (" 






dr 




f 


drr' 


xdm^{T 

a 


) / TT^ TT m 


a27r2 


JTc 




V12 4 r 







m 



Integrating and using just the first two terms inside the parenthesis, one gets 



(40) 



B = BnT-^ 



+ 



' 24^ 
vV 



2 2 

m T 



+ 



12a7r 



T 

Tc 



1 - 2 ) 

a 



drr am 



+ 



(41) 



On the other hand, for <^ 1, Eq. (!37|) is reduced to 

T V2 2V2n ^ ' 

At very high temperatures, the coupling constant, Eq. (1381) . decreases very slowly, implying 
that the ratio m/T is a slowly decreasing function at T/Tc ^ 1. Hence, m/T can be regarded 
as a constant, and the integral can be easily done. Rewriting the integrand in m/T powers, one 
can compute the integral of the B term at very high temperatures: 



B ^ BqT^ - 
= BoTi- 



vV 
24^ 
vV 



7 \ 7 T^- 
2 - - T- 



a/ T2 4 
uV 1 



1 

a J 



+ 



Uair 

3rr. 7 



12a 4 — 7/0; 4a7r4 — 7/0; a 

Note that Eq. fH51) is not valid for ^ = 4. For - = 4, one has 



7\ ^m^T^-^^ 
a) "t3 4- 



1 

a J 



(43) 



B ^ BoT^ 



[1 + 2\ogT)m^T^ + ^^{1 + 3\ogT)m^T, 



24a 



12a7r 



1 = 4. 

a 



(44) 
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As the asymptotic mass is proportional to T in first order, one has to maintain the Bq term. 
To obtain the asymptotic thermodynamics functions, we write the asymptotic B in terms of 
ideal massless gas pressure Pq: 



p /906o 15 n+21ogr)"^% ( 



p /, 90 a_4 15 1 
^0 % — -T-- - 



m 



+ 



45 



31ogT) 
1 



m 

3\ 



7 
a 



4; 



y z/7r^ 2Q;7r^ 4 — 7/0; ' 20:7?^ 4 — 7/0; ^ ' a' ' 

where 60 = The thermodynamics functions can now be computed. For instance, with the 
help of Eqs. dS]), dSI]) and fll^ one is able to obtain the asymptotic expression for the pressure 
for a 7^ and 7/0; = 4: 



P'(T) 



906o« , 15 

+ logT 



45 nr 



logT+... 



From Eq. ( l38i) . it is seen that the T logarithm is inversely proportional to g^^. Thus, with the 
help of the asymptotic mass Eq. (135|) . one obtains the asymptotic pressure: 



P\T) = Po 



^ 906o« 30 90 Qs 
1 — + — + 1=- + ... 



z^vr2 11 11^ -K 

When the coupling constant goes to zero, P'(T) should be that of an ideal massless gas. To this 
aim, one has to choose 60" = z^7r^/33. The pressure is then 



P'{T) = Po 



90 

1 + + 



2 = 4. 

a 



(46) 



The pressure for 7/0; 7^ 4 is easier to work with. Using Eqs. ([5]), (13T1) and (HSl) one gets 



P'iT) = Po 



15 2-^ /m 



7^ + 



15 1-? 



7^7 - "^0 



90 



T^-^ + ... 



2^4. (47) 
a 



47r24-^Vry 27r34-^VT/ z/tt 

a a 

In section II, the general mathematical solution for a consistent thermodynamics of massive 
gluons was handled. We now want to address the problem from the physics point of view: which 
solutions are physically relevant? 
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The first feature that one has to keep in mind is that the gluon plasma must behave as 
an ideal massless gas in the limit T —>■ oo. This property implies that the pressure has to be 
proportional to in this limit. As it was seen, the solution for a = and 7/a = 4 does not 
have a problem in the T — 00 limit. For the other solutions, one has to go back to Eq. (147!) 
and analyze the pressure. As m/T is proportional to powers of Qs, the terms in m/T do not 
disagree with the ideal massless gas limit. On the other hand, the term involving boa does. The 
only way to solve this problem is to make b^a term vanish, what implies that the powers of T 
must be negative. Hence, the general solution does not make physical sense for 7/a > 4. As 
a result, the general solution obtained in section II has possible physical meaning only if the 
condition 7/a < 4 is satisfied. We will call this condition of weak physical condition. Notice 
that solutions 1, 2 and 3 satisfy the weak physical condition. 

4 Analyzing the Physics Solutions 

We start with solution 1 at very high temperatures. The pressure in this scenario will be denoted 
by Pi, and can be written with the help of Eq. fHTl) and the condition 7 = as: 



Pi(T)=Po 



15 /m\^ 15 / 
1-^(7^ +7r^ 7^] +... 



^7r2 VT 



3 \TJ 



(48) 



Introducing the corrected asymptotic mass, Eq. fj42|) . in Eq. psj) . one obtains 



Pi = Po 



15 



^9s + 



15 



^9s 



'3^2 + 1^ 



+ 



(49) 



167r2^" ' 47r3^^ V 4v^ 
On the other hand, the perturbative QCD expression for the pressure at high temperature [21], 
known so far up to order g^\og{l/ gs)^, is: 



QCD 



16 TT^ 4 7T^ 



(50) 



Comparing Ens. (jUj) and fl50|) . one sees that the zero and second order terms match, but the 
term does noljj. Even so, besides the weak physical condition being satisfied, and the matching 



^The number inside the parenthesis which multiplies in the Eq. is ~ 0.93. 
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between asymptotic pressure and perturbative QCD up to g'^ order, solution 1 has another 
quahty: the connection between thermodynamics and statistical mechanics can be written in 
such a way which preserve the same form of standard statistical mechanics, Eq. f[T^ with 7 = 0. 

One has to remember that the extra term B is usually interpreted as a modification in the 
zero point energy, in other words, the vacuum is modified. If the vacuum changes, and now it 
is T dependent, it is expected an entropy associated with this vacuum, and consequently a T 
dependent vacuum entropy. But solution 1 does not have an extra term associated to a vacuum 
entropy, which can be seen as a limitation of this particular solution. 

The pressure for solution 2 at extremely high temperatures, denoted by P2, is: 



P2 =Pn 



1 - 



15 /m\ 



127r2 



+0(m7T^) + ... 



Notice that there is not a correction in m? /T^. This happens because the extra term B cancels 
out the term coming from the SM pressure. If one uses the HTL asymptotic mass in the pressure 
P2, then 



P9 = Pn 



1 - 



15^ 15^ 
24 7r2 24 vr 



+ . . . 



(51) 



The pressure P2 does not match with pQCD at any order in gg. The term g1 is about 1.5 
smaller than pQCD, while the gl is 6 times smaller. However, there is one way for solution 2 to 
match the pQCD result. If one uses m? = Sg'^T'^/A for the mass, instead of the HTL asymptotic 
mass, the perturbative result is recovered. In Ref . [12] , the gluon mass was taken as the plasma 



frequency 



UJ 



"^^"li^ • "^^^ effect of using the plasma frequency leaves the pressure 



more convergent than P2 and Pqcd- 

As it was seen in section 3.1, solution 3 is the only solution which has algebraic expressions 
for the thermodynamics functions. This solution, asymptotically, has a pressure P3 given by: 



Ps = Pn 



(52) 
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Comparing the asymptotic behavior of the P3, calculated with the HTL asymptotic mass, 
with the perturbative QCD pressure, one sees that there is no match at any order, as expected. 
However, if instead of the HTL mass one uses m = ^ for the mass, the QCD pressure up to 
is recovered. Motivated by the HTL gap equation, one can modify the mass relation in general 
and try a matching also at higher orders. One way to do this is dividing Eq. fl37j) by two, and 
redefining ^ — m. At lowest order the mass will be gsT/2. Extending the mass equation to 
second order, the term for the pressure will match the perturbative QCD result well. Hence, 
as a bonus, solution 3 with a modified gap equation for the thermal mass, results in a match 
between P3 and Pqcd up to the g^ order. 

As the internal energy and entropy are changed by a T dependent function in solution 3, 
the extra B term can be regarded as a zero-point energy, with an entropy associated with it. 
It is important to emphasize that as one has already obtained an algebraic expression for the 
pressure at the whole temperature range, it is possible to fit m(T) by QCD lattice data, and 
search for deviations of the HTL asymptotic mass near of Tc. 

The solution 7/0 = 4 has also to be analyzed. This solution is the only one which has a 
term of order gg. Although it satisfies the weak physical condition, the contribution of order 
gs contradicts any pressure computation in the literature, either HTL, 2-Loop ^-derivable 
approximation or perturbative QCD. The absence of a term in gg can be justified by thermal 
field theory, since the first-order correction to the partition function is proportional to g^ [26j. 
Therefore, we will let this solution out of the possible physical solutions. 

5 Discussion and Summary 

In this work we have studied the thermodynamics and the statistical mechanics of a system with 
a temperature dependent Hamiltonian in the canonical ensemble. Earlier works have shown that 
a system with a T dependent Hamiltonian is thermodynamically inconsistent if the connection 
between thermodynamics and statistical mechanics is the same as in the T independent Hamil- 
tonian case. In order to have a general solution to this problem, we developed a formalism which 
gives a general connection between thermodynamics and statistical mechanics, where all the 
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thermodynamics relations are satisfied. We have seen that it is necessary to add an extra term 
to the thermodynamics relations to guarantee the thermodynamics consistency, and depending 
on the choice that is made for the parameters a, 7 and 1], the extra term B can be simpler. 

The general solution developed here is a generalization of the solution proposed by Goren- 
stein and Yang [17]. In their work, they have introduced one particular manner to maintain 
the thermodynamics consistency: a correction coming from the zero point energy of the Hamil- 
tonian modifies the pressure and the internal energy of the system, while leaving the entropy 
unchanged. In the present work, it is shown that letting the Hamiltonian function to be depen- 
dent on the temperature of the system, represented by a function /(T), implies that there are, 
in principle, a large number of ways to render the thermodynamics of the system consistent. 
We emphasize that this procedure is applicable in any system with a T dependent Hamiltonian, 
and not only in the case of quasi-particle models. 

For an ideal quasi-gluon gas, we computed the pressure, the energy and the entropy, within 
the formalism developed in section 2. This calculation, named solution 3, is easier to work with 
as no integral in T is necessary to determine the extra term B. As the formalism does not provide 
an expression for the thermal mass m{T), we used the mass calculated in the HTL approach. 
We achieved an optimal fit for the thermodynamics functions with two free parameters in the 
expression for the strong coupling at finite T. 

We also studied the general asymptotic behavior of the solutions contained in a 7^ cases. 
The analysis of these solutions resulted in a reduced number of possible physical solutions. The 
main properties of solutions 1, 2, and 3 were discussed, and it was showed that solution 1 is 
the only one of the three solutions which recovers pQCD to order g"^ when the HTL asymptotic 
mass is used. In order for solutions 2 and 3 recover pQCD, it is necessary to change the gap 
equation for the mass. Finally, we saw that the only solutions that have physics meaning are 
those where the constants a, 7 and rj satisfy the weak physical condition, a = or ^ /a < A, 
together with the constraint a = 7 + r^. An extension of the present calculation to the case of 
a finite chemical potential will be presented in the near future. 

This work was supported by FAPESP (04/15276-2) and CNPq (307284/2006-9). 
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A Appendix 



We need results for integrals like 

, , 1 /■°° , x""^ 1 , , 

In[r) = ^ , , / dx 1 1 . (53) 

T{n) Jo (x2 + r2)l e(^'+'-')^ - 1 

The relevant integrals for our calculation, following [211 EZj and generalizing their results, are 
the following: 

7]-2 ~- / r 1 \ 1 °° 

^3(-) = - - log + 7i. - d - E -mim + 2)r^('"+i). (54) 



223 22 23V 227r 2 J 2\^^^ 

hir) = ^,-r'f, + f/ + "^ (log + - I) + ^ E (55) 

where = 23m+i"2"m('^^^^|^"^'' • If one tests the series for convergence, for instance the ratio 
test, one finds that the series are convergent for r < 27r. 
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